Quantum Fluctuations and Excitations in Antiferromagnetic Quasicrystals 
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We study the effects of quantum fluctuations and the excitation spectrum for the antiferromag- 
netic Heisenberg model on a two-dimensional quasicrystal, by numerically solving linear spin-wave 
theory on finite approximants of the octagonal tiling. Previous quantum Monte Carlo results for the 
distribution of local staggered magnetic moments and the static spin structure factor are reproduced 
well within this approximate scheme. Furthermore, the magnetic excitation spectrum consists of 
magnon-like low-energy modes, as well as dispersionless high-energy states of multifractal nature. 
The dynamical spin structure factor, accessible to inelastic neutron scattering, exhibits linear-soft 
modes at low energies, self-similar structures with bifurcations emerging at intermediate energies, 
and flat bands in high-energy regions. We find that the distribution of local staggered moments 
stemming from the inhomogeneity of the quasiperiodic structure leads to a characteristic energy 
spread in the local dynamical spin susceptibility, implying distinct nuclear magnetic resonance spec- 
tra, specific for different local environments. 

PACS numbers: 03.75.Hh,03.75.Lm,05.30.Jp 



I. INTRODUCTION 

Quantum fluctuations are responsible for various de- 
grees of disorder in low-dimensional quantum antiferro- 
magnets. In particular two-dimensional systems show a 
variety of quantum disordered phases, competing with 
conventional long-range magnetic order. For example, 
while the Heisenberg model on the square lattice exhibits 
true long-range order at zero temperature!., spatial inho- 
mogeneous magnetic exchange eventually leads to a com- 
plete suppression of magnetic order in structures such as 
the plaquette lattice, driven by local singlet formation^. 
Other sources of magnetic disorder are frustration ef- 
fects due to competing interactions*, and site/bond de- 
pletion^, where the reduction in the long-range magnetic 
order is accomplished by proliferation of localized low- 
energy excitations 6 . All these systems share the trans- 
lational invariance of the underlying lattice structure - 
assuming sufficient self averaging in the case of quenched 
disorder. 

Quasiperiodic systems, lacking translational symmetry 
in addition to their inhomogeneous lattice structure were 
initially thought to not support sizeable correlations of 
localized magnetic moments. However, recent neutron 
scattering experiments on Zn-Hg-Hc 7 as well as Cd-Mg- 
Tb* icosahedral quasicrystals support the presence of sig- 
nificant magnetic correlations in these three-dimensional 
quasicrystalline compounds, with the absence of true 
long-range magnetic order— due to large frustrations in 
the antiferromagnetic exchange 10 . 

Bipartite, and thus unfrustrated, quasiperiodic crys- 
tal structures were indeed shown to allow for sizeable 
two-sublattice antiferromagnetic order in a recent quan- 
tum Monte Carlo study of the spin- 1/2 Heisenberg model 
on the octagonal tilingi°«. Furthermore, the magnetic 
order in this two-dimensional system was found to ex- 
hibit nontrivial patterns in the local staggered moment 




FIG. 1: Finite approximant of the octagonal tiling with 239 
sites (thin lines), along with a superimposed inflated 41 sites 
approximant (thick lines) of rescaled edge length by a factor 
of X = 1 + v2- For the 41 sites approximant the vertices of 
one of the two sublattices are dressed with a disk to exhibit 
the bipartite nature of the octagonal tiling. 



distribution, reflecting the self-similarity of the underly- 
ing quasiperiodic lattice structure^ -. A renormalization 
group approach based on this self-similarity indeed gives 
a gross account on the observed spread in the staggered 
magnetization^*. Furthermore, the static spin structure 
factor exhibits magnetic selection rules that impose a 
shift of reciprocal space indices* -. The resulting neu- 
tron diffraction pattern* - can be accounted for by anal- 
ysis^ of the quasicrystal spin groupie, also applicable to 
frustrated classical models on the octagonal tiling with 
long-range exchange interactions**. 

While static properties of the Heisenberg antiferromag- 
net on the octagonal tiling are thus well studied, little is 
known about the spectral properties of these systems. 



2 



On general grounds one would expect gapless Goldstone- 
modes to dominate at low-energies, even though the 
translational symmetry is absent. In addition, one would 
expect to find multifractional eigenstates, as observed in 
tight-binding models on quasiperiodic lattices^. Here, 
we investigate dynamical properties of quantum mag- 
netic quasicrystals in order to identify the relevant energy 
scales of quantum fluctuations in such systems, determine 
the magnetic excitation spectrum, as well as to provide 
theoretical grounds for future experiments on magnetic 
quasicrystals, such as inelastic neutron scattering or mag- 
netic resonance. In particular, we use linear spin-wave 
theory, which was successfully used in studies of peri- 
odic magnetically ordered systems, and apply it to the 
quasiperiodic case. 

The outline of the paper is as follows: Basic properties 
of the octagonal tiling are presented in the following sec- 
tion, in sec. nni we review linear spin wave theory in a 

real space formulation, and present a numerical construc- 
tion of the eigenmode expansion. An alternative scheme 
is given in the Appendix. The results of applying this 
method to the octagonal tiling are discussed in Sec. IIVI 
In Sec. lIV Al we discuss static properties of the magnetic 
correlations on the octagonal tiling, and compare our re- 
sults to previous quantum Monte Carlo simulations. A 
detailed analysis of the excitation spectrum is presented 
in Sec. IIVBI followed by a discussion of dynamical mag- 
netic properties, such as the dynamical spin structure 
factor in Sec. II V CI and the local dynamical spin suscep- 
tibility (Sec. IIVD|l . Finally, we conclude in Sec. IVl with 
a perspective on future investigations. 

II. OCTAGONAL TILING 

In the following, we analyze the magnetic ground 
state properties and excitation spectrum of the nearest- 
neighbor antifcrromagnetic spin— 1/2 Heisenberg model, 

H = jJ2SiSj, J>0, (1) 

(hi) 

on the most prominent example of a magnetic quasicrys- 
tal in two dimensions, the octagonal tiling. The octago- 
nal tiling is a bipartite quasiperiodic crystal system, and 
possesses an overall eightfold rotational symmetry, allow- 
ing for simple two-sublattice antiferromagnetismi^. Sites 
in this tiling have coordination numbers z ranging from 
3 to 8, leading to a broad distribution of local staggered 
moments in the magnetically ordered ground stated. A 
further important property of the octagonal tiling, in the 
absence of translational invariance, is its self-similarity 
under inflation transformations^. This reversible oper- 
ation refers to a well-defined decimation of a subset of 
vertices of the tiling, followed by a re-connection of the 
new vertices. Aside from a trivial rescaling of the length 
scale by a factor A = 1 + the infinite quasicrystal is 
left unchanged by this transformation. 



For our numerical study we consider finite square ap- 
proximants of the octagonal tiling with 41, 239, and 1393 
sites. These approximants can be obtained by the "cut- 
and-project" method from a four-dimensional cubic lat- 
tice^, and are related by the inflation transformation, an 
example of which is shown in Fig. ^ In order to avoid 
boundary-induced frustration effects, we apply toroidal 
boundary conditions^. Due to the lack of translational 
symmetry, we thus need to solve real-space linear spin- 
wave theory on lattices with up to 5572 sites. Before 
presenting our results, we provide details about the nu- 
merical scheme used in our calculations in the following 
section. 



III. NUMERICAL SPIN- WAVE 
APPROXIMATION 

In this section, we review linear spin-wave theory, ap- 
plied to the antiferromagnetic Heisenberg model on finite, 
bipartite lattices. We consider a bipartite lattice with 
sublattices A and B consisting of Na and Nb sites, re- 
spectively. For the octagonal approximants considered in 
this work Na = Nb , however the following approach also 
applies if Na and Nb are different. Following the stan- 
dard Holstein-Primakoff approach^, we represent the 
spins in terms of bosonic operators and bi. For i G A, 

S* = S - atoi, 

_/ t \ 1/2 
S+ = V2^(l-^J a h (2) 

/ t \ 1/2 

*--V5*}(i-as) , 

where S denotes the spin magnitude, and 

s; = -s+bfa, 

/ t \ 1 I 2 
St = V2Sb] (l-^) , (3) 

for j 6 B. The linear spin- wave Hamiltonian is obtained 
by substituting the above identities in Eq. |JT]|. expanding 
the square roots in of 1/5, and keeping terms of lowest 
order (1/S)°, 

Hsw = -JS(S+l)N b + JSH 2 , (4) 
H2 = J2 (ajai + tyfej + atftt + & jai ) . ( 5 ) 

where Nb denotes the number of bonds, and the sum in 
i?2 extends over all bonds of the bipartite lattice. Intro- 
ducing the N s = Na + Nb component row vector 

a) = (a\, ...,a\ !A ,b 1 , ...,b NB ), (6) 
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and the corresponding columnar conjugate, we can ex- 
press the quadratic part of the Hamiltonian in matrix 
notation, 

H 2 = a) Ma = ^ a^M^a, (7) 

where 

(m^A = 6 l , k 5 i .i+6 j+NA . k S J+NA ,i 

V / k,l 

+5i t k$j+NA,l + fij+NA,k$i,l (8) 

is the connectivity matrix of the lattice structure. The 
bipartiteness of the lattice thus allows for a direct for- 
mulation in terms of the N s x N s hermitian matrix M, 
instead of the general formulation based on a 2N S x 2N S 
matri»2i. We now seek a Bogoliubov transformation to 
the N n normal bosonic modes [3k, such that the spin- 
wave Hamiltonian Hsw is diagonal when expressed in 
terms of the @k, 

H sw = N s E + J2 u> kPtPk, (9) 
fc=i 

where Uk > denotes the eigenfrcqucncy of the fc-th 

mode, and E the ground state energy per site. To this 
end we make an Ansatz, 

a = T0, (10) 

with 

P = Wt,...,0 t N+ ,0N + +l,:.,0Nn), (11) 

so that we divide the N n normal bosonic modes into two 
disjunct sets of length N + , and iV_ = N n — N + , re- 
spectively, further determined below. Due to the bosonic 
commutation relations, the transformation T has to fulfill 

rrrt = e, (12) 

with matrices T and E, defined by 




Here, ljv denotes the N x N identity matrix-''. Since 
E 2 = 1, as well as T 2 = 1, and because the left and right 
inverse of a square matrix are identical, if T satisfies 

T f ET = T, (14) 

the bosonic commutation relations are fulfilled. Since T 
diagonalizes H2, we have 

T^MT = fl, fi = diag(wi,...,WArJ, (15) 

from which we find upon multiplication from the left with 
IT, that T has to satisfy 

EMT = TTO. (16) 



The column vectors of T are thus seen to be related to the 
right eigenvectors of EM. For semi-positive M, EM has 
real eigenvalues, and eigenvectors belonging to different 
eigenvalues are orthogonal 2 ^. We denote the number of 
positive (negative) eigenvalues by N + (N-), the number 
of zero- modes by No, and label the positive (negative) 
eigenvalues by Xf , i — 1,...,N + (A7, i — l,...,iV_). Af- 
ter numerically solving the non-hermitian N s x N s eigen- 
value problem for EM 22 , we construct within the sub- 
space of each degenerate eigenvalue > 0(< 0) of di- 
mension eigenvectors z^ 1 , d ± , obeying 

&^nj = ±*J, l,J = l,-,dt (17) 

using Gram-Schmidt orthogonalization with respect to 
E (ref. 123) . We then obtain T from the orthonormal 
eigenvectors as 

T = «i, 4 +>d + + , z~ di ,..., z-_ >d __ ), (18) 

where Eq. l(T7|) ensures that Eq. lfTl|l is satisfied. The 
corresponding eigenfrequencies are given by 

uj k = JSX+, k = l,...,N+, (19) 
uj N+ +k = -JSX k , fc = l,...,JV_. (20) 

The zero-modes of EAf correspond to collective modes 
due to the broken continuous symmetry implied by the 
classical Neel stat o 6 i 21 . Furthermore, the ground state 
energy becomes 

N- 

E N s = -JS(S+l)N b + JSj2\ X k\- ( 21 ) 

fe=i 

An alternative means of numerically constructing the 
transformation matrix T is presented in the appendix. 
We have verified, that both approaches indeed yield the 
same results. 



IV. RESULTS 

In this section, we present the results obtained by ap- 
plying the linear spin-wave approximation to the octag- 
onal tiling introduced in Sec.^ Within linear spin- wave 
theory, ground state expectation values of both static and 
dynamic magnetic correlations can be calculated from 
contractions of bosonic normal mode operators. It turns 
out convenient for this purpose, to express the Bogoli- 
ubov transformation in terms of the row vectors 

a 1 = (a[,...,a , NA ,b[,...,V NB ,a 1 ,... 1 a NA ,b 1 ,...,b NB ), 

p = (/?!,..., 0* Nn , 0i,.. .,p Nn ), 

as 

a = fp, f = [!L M . (22) 
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One obtains U and V from the transformation matrix T 
of Eq. H1U|) . upon denning N + x N + matrices A and B, 
and iV_ x N- matrices C and D such that 

r = (i g) , (23) 

and gets 

(o D °-) ■ ^ S) • <-) 

A. Staggered magnetization 

In this section we examine static properties of the mag- 
netic ground state in the octagonal tiling. For this pur- 
pose, we first calculate the staggered magnetization at 
each lattice site i in linear spin- wave theory, 

m s (zH|(Sf)| = S-£|^ fc | 2 . (25) 

k 

The spatially averaged staggered magnetization 

m a = m s {i) (26) 

i=l 

is shown as a function of the system size N s of the approx- 
imant of the octagonal tiling in Fig. El The finite size val- 

— 1/2 

ues scale well as a function of N s , and indicate a size- 
able staggered magnetization of the long ranged ordered 
ground state in the octagonal tiling. The value of m s 
extrapolated to the thermodynamic limit is m s = 0.34, 
and agrees well with the quantum Monte Carlo result, 
m s = 0.337 ± 0.003iS. The weak reduction of the or- 
der parameter by quantum fluctuations is indicative for 
the feasibility of the linear spin-wave approach in this 
system. The ground state energy Eq, shown in the in- 
set of Fig. [21 is also found to scale well as a function of 
N s 3 ^ 2 , as expected for an ordered stated. The extrapo- 
lated value in the thermodynamic limit, Eq = —0.646, 
compares well with the quantum Monte Carlo result, 
E Q = -0.6581(1)2^. 

In Ref. 0, the magnetic ground state on the octago- 
nal tiling was found to exhibit a nontrivial local structure 
reflecting the self-similarity of the underlying quasiperi- 
odic lattice structure. We now analyze to what extent 
spin-wave theory is able to reproduce this structure and 
thus to account for the specific nature of quantum fluc- 
tuations in an inhomogeneous connectivity. For this pur- 
pose, we show in Fig.|3]the linear spin-wave results of the 
local staggered magnetization, m s (i), for the 1393 sites 
approximant, and compare those with values grouped ac- 
cording to the coordination number z of the various sites, 
with the quantum Monte Carlo data of Ref. fiot 

We find linear spin- wave theory to qualitatively repro- 
duce characteristic features of the local staggered mo- 
ment distribution, such as (i) a wide spread of the mo- 
ments, in particular for small values of z, (ii) a prominent 
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FIG. 2: Finite size scaling of the ground state staggered 
magnetization, m s , and the lowest excitation gap, A, for the 
5 = 1/2 Heisenberg model on the octagonal tiling in linear 
spin-wave theory. The inset shows the finite size scaling of 
the ground state energy, Eq, for the same system. 



bimodal splitting of the moments for z = 5 and (iii) the 
hierarchical structure observed in the splitting of the mo- 
ments for sites with z = 8, shown in the inset of Fig. |3] 
These splittings in the local staggered moments can be 
accounted for by the properties of inequivalent sites un- 
der deflation transformations, reflecting their different 
local environments^. In particular, fivefold sites with 
z = 5 always occur in pairs, with two different types of 
site. The first type is connected to four fourfold (z — 4) 
sites and the other fivefold site, while for the other type 
two neighbors are fourfold, two threefold (z = 3) and one 
fivefold. This difference in the local connectivity leads to 
the observed splittingiS, and is also reflected in a different 
behavior of the two types of fivefold sites under deflation 
transformations: while one type is decimated, the other 
remains as a threefold site in the deflated tiling. The 
eightfold (z = 8) sites exhibit an even richer, hierarchi- 
cal structure shown in the inset of Fig. |3J with moments 
grouped according to the different deflation properties of 
the eightfold sites, namely their new coordination num- 
ber z'. The high symmetry sites with z' = 8 show a fur- 
ther hyperfine splitting, with moments grouped accord- 
ing to the value of z" under a second deflation, indicated 
by numbers next to the symbols in the inset of Fig. [3] 
As discussed in Ref. 0, these splittings eventually lead 
to a multifractal distribution of local staggered moments 
in the infinite quasicrystal, for this class of sites. 

From a quantitative comparison of the local staggered 
magnetization between linear spin- wave theory and quan- 
tum Monte Carlo, we find characteristic limitations of the 
spin-wave approach to persist on a local level. Namely, 
while local staggered moments of high-connectivity sites 
are reproduced even quantitatively, deviations of about 
8% are observed for low-connectivity sites. In an inho- 
mogeneous environment, linear spin- wave theory is thus 
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FIG. 3: Dependence of the local staggered magnetization 
on the coordination number z for all sites in the 1393 sites 
approximant of the octagonal tiling within linear spin-wave 
theory of the S = 1/2 Heisenberg model (LSWT). For ref- 
erence, results of quantum Monte Carlo simulations of the 
same system (QMC) are also shown—. The inset exhibits the 
hierarchy of the local staggered magnetization of the z — 8 
sites, grouped according to the value of z' under a deflation 
transformation. Numbers next to symbols give the value of 
z" for z' = 8 sites under a further deflation. 



more accurate at sites of large coordination, as might 
have been expected from its behavior in homogenous sys- 
tems. 

The long-range antiferromagnetic order in the octago- 
nal tiling leads to characteristic neutron diffraction pat- 
terns, due to selection rules imposed by the magnetic 
symmetrjii2ii£. These patterns can be obtained from the 
static longitudinal structure factor 



S'l(k) 



zk- (r$— t» ) 



(27) 



which within linear spin-wave theory amounts to the 
Fourier transform of the real space distribution of (S*). 
Since linear spin- wave theory reproduces the spatial stag- 
gered moment distribution rather well, as seen from 
Fig. |3 the resulting static longitudinal structure factor, 
shown in the left of Fig.0]also compares well to the quan- 
tum Monte Carlo result. In particular, it exhibits the 
extinction of nuclear Bragg peaks and the emergence of 
new, magnetic Bragg peaks. For a detailed theoretical 
derivation of the various selection rules in the octa gon al 
tiling, we refer to the explicit enumerations in Ref. Il2l 

The structural distribution of local staggered moments 
observed in Fig.[3]exhibits an inhomogeneous distribution 
of quantum fluctuations. As seen from Eq. I|25l) . fluctua- 
tions in m s arise from distinct contributions from the var- 
ious eigenstates of the system. To quantify the relevance 
of the different eigenstates for the quantum fluctuations, 
we study the reduction of the staggered magnetization 
at each lattice site i, which can be parameterized as a 
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FIG. 4: Intensity plot of the static longitudinal magnetic 
structure factor ^'(k) (left), and the integrated dynamical 
spin structure factor S x (k) (right) for the S = 1/2 Heisen- 
berg antiferromagnet on the 1393 sites approximant of the 
octagonal tiling. 



function of energy, 



J2 k S(U - LOk) 



(28) 



In Fig. [5] we show Sm s (i,u>), averaged separately over 
sites with z — 3, and z = 

In both cases, the dominant contribution to quan- 
tum fluctuations stem from the low-energy modes, with 
uj/JS < 2. For threefold sites, further contributions to 
8m s arise from higher energy modes, which are not rel- 
evant for the eightfold sites. In general, we find the up- 
per bound of the energy range that is relevant for quan- 
tum fluctuations to decrease with increasing coordination 
number z. Although sites with low coordination numbers 
thus receive quantum fluctuations over a larger range in 
energy space, their staggered moment is typically larger 
than for high coordinated sites, as seen from Fig. [21 A 
transfer of relevant quantum fluctuations to low-energy 
modes is thus responsible for a decrease of the staggered 
moment at specific sites. This increased relevance of low- 
energy modes for quantum fluctuations is also observed 
for sites with the same coordination number. For exam- 
ple, the inset of Fig. shows 6m s for the fivefold sites, 
where the bimodal splitting was observed in Fig. [3J av- 
eraged separately over sites with a small and a large mo- 
ment, respectively. The additional reduction of the local 
moment for one type of fivefold sites is clearly seen to 
be due to a proliferation of quantum fluctuations in the 
lower energy region. 

This hints at a close link between the structure of the 
inhomogeneous magnetic ground state and the magnetic 
excitation spectrum. In the following subsection, we pro- 
ceed to analyze the magnetic excitations in more detail. 



B. Excitation Spectrum 

In the previous section we found that spin-wave theory 
provides a qualitative, and even quantitative account on 
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FIG. 5: Frequency dependence of quantum fluctuations to 
the staggered magnetization, 5m a {ui), for the 5=1/2 Heisen- 
berg model on the octagonal tiling, averaged separately over 
sites with coordination numbers z = 3 and 8. The inset shows 
dm s {w) for z = 5, averaged separately over sites with a small 
and large staggered magnetization, respectively. 
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FIG. 6: Excitation spectra of the Heisenberg antiferromagnet 
on finite approximants of the octagonal tiling of different sizes 
within linear spin-wave theory. The inset shows the density 
of state, p(io), for the largest approximant. 



the magnetic ground state properties of the Heisenberg 
model on the octagonal tiling. We now proceed to use 
linear spin-wave theory to gain insights into the spectral 
properties of this quasiperiodic antiferromagnet. 

The spectra of spin-wave normal modes for the first 
three approximants of the octagonal tiling are shown in 
Fig. For finite approximants the spectra consists of dis- 
crete sets of energy levels, spanning an extended range 
up to uj max / JS w 7.3. This implies a bandwidth almost 
a factor two larger than for the Heisenberg model on the 
square lattice, where uj max / JS — 4. Upon increasing the 
system size, the energy spectrum appears to develop a 
dense band, and two isolated pockets at higher energies, 
near ui/JS ~ 6.5 and 7.3, respectively. The presence of 
a third pocket near uj/ JS ~ 5.5 cannot be excluded from 
the finite size data. However, from Fig. the gap near 
ioj JS w 5 appears to eventually close for higher approxi- 
mants, whereas the gaps to the higher energy pockets do 
not show any decrease for the approximants considered 
here. 

For a more quantitative analysis of the distribution of 
normal modes, we calculate the density of states of the 
spin- wave excitation spectrum (DOS), 

p{w)=^5(w-w k ), (29) 
k 

shown for the largest approximant in the inset of Fig. 
The spin- wave DOS exhibits a characteristic spiky shape 
in the high-energy regime, in particular for energies larger 
than u/JS » 3, similar to shapes found for the tight- 
binding Hamiltonian on the same lattice structure^!. In 
the low-energy region, the DOS appears more smooth, 
with a residual roughness due to the limited resolution 
of the energy spectra due to finite size effects. We indeed 



find similar resolution limited roughness also on finite 
square lattice systems. 

On general grounds, the presence of a long-range or- 
dered ground state is expected to characterise the low- 
energy properties of the excitation spectrum. In partic- 
ular, due to the broken SU(2) symmetry, we expect the 
gap to the lowest excitation, A, to close in the thermo- 
dynamic limit of the quasiperiodic tiling. A finite size 
scaling analysis of A, shown in Fig. [3 is indeed consis- 
tent with A cx N~ a , where a « 0.5. At low energies the 
DOS furthermore increases linearly, p(uS) cx to, as seen 
from the quadratic low-energy behavior of the cumula- 
tive density of states, 

N(u)= r<Me), (30) 

To 

shown in Fig- El Here, we employ N(u>), since this quan- 
tity is less susceptible to resolution limited roughness 
than p(u>), due to the frequency integration. The low- 
energy features of the spin-wave DOS in the octagonal 
tiling are also observed in periodic antiferromagnets. For 
example, on the square lattice, antiferromagnctic spin- 
waves obey a linear dispersion relation at low energies 
Wk rs v / 8^5 , |k|, resulting in a linear low-energy DOS in 
this system. 

The similarity of the low-energy DOS for the Heisen- 
berg model on the quasiperiodic octagonal tiling to the 
periodic case suggest that the low-energy excitation could 
be delocalized, magnon-like modes also in the octagonal 
tiling. We are thus lead to analyze the spatial extent 
of the eigenstates found within linear spin-wave theory. 
An appropriate method to characterize the localization 
properties is the inverse participation ratio (IPR) 28 which 
expresses the spatial extent of the wavefunction of a given 
state. For the bosonic excitations of quantum magnets, 
the IPR of an eigenstate is given by its contribution to 
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FIG. 7: Cumulative density of states, N(u), of the Heisen- 
berg antiferromagnet on finite approximants of the octagonal 
tiling of different sizes within linear spin-wave theory. The 
inset exhibits the quadratic scaling of N(lj) at low energies 
for the largest approximant, indicative of a linear low-energy 
density of states. 



the quantum fluctuations of the staggered magnetization. 
Following Ref. 0, we define the energy dependent inverse 
participation ratio 



J( W ) 



J2k !k Sjco - u k ) 



where 



h = 



(EJW) 2 



(31) 



(32) 



and study its scaling behavior upon increasing the sys- 
tem size. For delocalized states of a d dimensional quan- 
tum system, the IPR decreases with the system size, -/V s , 
as N~ . Exponentially localized states should be very 
insensitive to the system size and one expects a size in- 
dependent IPR. 

Fig. [8] shows the calculated IPR as a function of en- 
ergy for the spin— 1/2 Heisenberg model on finite approx- 
imants of the octagonal tiling. For comparison, we show 
in the inset the corresponding quantity for the Heisen- 
berg model on the square lattice. In the square lattice 
case, all eigenstates show the characteristic N~ x scaling, 
expected for extended magnon states. For the octago- 
nal tiling, we find such behavior only for the low-energy 
modes, but due to the limited resolution cannot exclude 
a reduced finite size scaling down to zero frequency. The 
low-energy modes in the octagonal tiling thus appear as 
extended excitations out of the antiferromagnetic ground 
state, similar to coherent magnons in the periodic case. 
The higher energy states also do not appear exponen- 
tial localized, but show a significantly reduced finite size 
scaling I(uj) cx N~@, with a scaling exponent < 1, 
characteristic for multifractal states, as observed for crit- 
ical states at the Anderson localization transition 29 . The 
effective exponent f3 decreases towards the upper edge of 
the spectrum, being lowest for energies near the isolated 
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FIG. 8: Frequency dependent inverse participation ratio, 
I(u)), for the S — 1/2 Heisenberg antiferromagnet on finite 
approximants of the octagonal tiling of different sizes, ob- 
tained within linear spin-wave theory. Dashed lines indicate 
a finite size scaling as iV" 1 , expected for extended states. The 
inset shows linear spin- wave results for on square lattices 
of different sizes. 



pockets found in the DOS in Fig. HJ1 While a reliable 
determination of (3 would require the study of substan- 
tially larger approximants, we estimate a value of /3 ~ 0.2 
for the high-energy modes. The multifractality of eigen- 
states found in non-interacting models on quasiperiodic 
crystals^ is thus clearly observed also for excitations of 
strongly correlated systems, such as the Heisenberg an- 
tiferromagnet. 



C. Dynamical spin structure factor 

Having analyzed the spin-wave excitations of the 
Heisenberg antiferromagnet on the octagonal tiling, we 
proceed to study magnetic properties exhibiting the mag- 
netic excitation spectrum, such as the dynamical spin 
structure factor, and local dynamical spin susceptibility. 

The dynamical spin structure factor, accessible ex- 
perimentally by inelastic neutron scattering, reflects the 
time-dependent spin-spin correlation functions, trans- 
formed to momentum space, 

1 Ns 

S x (k, W ) = — £ j^-^S^iJ,^, (33) 



with 



SHhj>") = 2 j_ *e**W(t)Sr(o) + Sr(t)S+(0)). 

(34) 

Using the normal-mode expansion, Eq. Q22[l. we obtain 
the following expression for S (i, j, u>) within linear spin- 
wave theory, 

= S^ (U ik U* k + V* k V jk ) 6{u - Wfc ), (35) 
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FIG. 9: Intensity plot of the dynamical spin structure factor, 
£(k,w), for fixed lo/JS = (left) and w/JS = 1.8 (right) 
for the 5=1/2 Heisenberg antiferromagnet on the 1393 sites 
approximant of the octagonal tiling. 



from which S (k, oj) is readily obtained using a fast 
Fourier transformation. Most spectral weight in <S lJ "(k, ui) 
is located at momenta corresponding to magnetic Bragg 
peaks, as seen by comparing the integrated dynamical 
spin structure factor, 



S x (k) = /^* x CM0 



(36) 



shown in the right of Fig. 0] to the static longitudinal 
structure factor, 5"(k), in the left of Fig. For a de- 
tailed analysis of the dynamical spin structure factor, we 
consider both constant-frequency scans, shown in Figs. [21 
to 1121 as well as scans along various momentum space 
directions, shown in Figs. El an d El To increase the 
contrast in these figures, we have rescaled the data to 
the maximum value in each plot, separately. 

In order to compare the relative spectral weight at 
different energies, we show in Fig. 1101 the momentum- 
integrated spectral function, 



d 2 k 



S(k,Lj). 



(37) 



Compared to the DOS, we find that apart from the low- 
energy region below oj/JS s» 2, the shape of <S , - L (w) 
closely reflects the structures in the DOS. In the low- 
energy region, we observe an disproportionately large 
contribution to S (to), in contrast to the low DOS in 
this region. The difference between ^(w), and the DOS, 
p(lo), in linear spin- wave theory is obtained using Eq. 
(EH as 



S ± (u) - Sp(co) = 2SJ2 \V ik \ 2 S(u ~ 



(38) 



and indicates, that the extra contributions to S (w) at 
energies u> / JS < 2 are due to the spectral predominance 
of low energy quantum fluctuations, 5m s {us), Eq. (|28|l . 
shown in Fig. In the following, we first consider the 
dynamical spin structure factor in this low-energy region, 
which is dominated by magnetic Bragg scattering. 




FIG. 10: Frequency dependence of the momentum-integrated 
dynamical spin structure factor, S ± (lo), for the S = 1/2 
Heisenberg antiferromagnet on the octagonal tiling. 
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FIG. 11: Intensity plot of the dynamical spin structure factor, 
S(k,w), for fixed lo/JS = 2.9 (left) and lo/JS = 3.0 (right) 
for the 5=1/2 Heisenberg antiferromagnet on the 1393 sites 
approximant of the octagonal tiling. 



In the elastic limit of S f± (k, lo — ► 0), shown in the left 
of Fig. 1^1 we can indeed identify the magnetic Bragg 
peak positions of the static longitudinal structure fac- 
tor, shown in Fig. 01 Furthermore, for all energies up to 
lo/JS w 2, similar patterns as in 5 _L (k, 0) are observed, 
albeit with the width of the peaks increasing upon in- 
creasing the energy. Eventually, these peaks evolve into 
ring-like structures, centered around the magnetic Bragg 
peaks, such as shown for lo/JS = 1.8 in the right part of 
Fig- El This is a clear indication for magnetic soft-modes 
at low energies, which dominate the magnetic response 
in this energy regime. 

Typical examples of 5 ,J "(k) at higher frequencies are 
shown in Fig. El and 1 121 We find all plots to exhibit an 
eightfold overall symmetry, as expected for the octago- 
nal tiling. However, the positions of the dominant peaks 
are different from the magnetic Bragg peaks found be- 
low lj/JS w 2. For example, for lo/JS = 3.0 (right of 
Fig- EJ) most spectral weight is located at k « (0, 1.37r), 
and k w (0, 4.37r), as well as symmetry-related momenta. 
Peaks at these momenta are absent in both S 1 " (k, 0), and 



FIG. 12: Intensity plot ol the dynamical spin structure factor, 
£(k,w), for fixed uj/JS = 4.6 (left) and lu/JS = 6.4 (right) 
for the 5=1/2 Heisenberg antiferromagnet on the 1393 sites 
approximant of the octagonal tiling. 



the nuclear Bragg scattering (Fig. 5 (a) of Ref. ITol). In 
fact, we do not observe pronounced spectral weight at 
nuclear Bragg peak positions for any finite energy cut: 
to give a further example of this fact, we show 5 _L (k) at 
uj/ JS = 4.6 in the left of Fig. Upon changing the 
energy-level of the cut only slightly, the patterns found 
in the high-energy region change more drastically than 
those in the low-energy regime. As an example, in the 
left of Fig.El we show 5 J -(k) for uj/JS = 2.9. Compared 
to S x (k) at uj/JS = 3.0 (right of Fig.HT}, we indeed find 
a different set of dominant peaks. 

The patterns in S l_L (k) for energies corresponding to 
the high-energy pockets in the DOS (c.f. Fig. [fjj) con- 
sist of more diffusive structures than those at lower en- 
ergies. For example, we find broad ring-like structures, 
centered around the magnetic Bragg peaks in S (k) for 
uj/JS = 6.4, shown in the right of Fig. The ap- 
pearance of such diffusive structures at high energies is 
expected from the results of Sec. IIVBI where the the 
high-energy excitations were found to be spatially less 
extended than the low-energy states. This reduced spa- 
tial extent leads to smeared diffraction patterns observed 
in Fig. [H 

In order to analyze the momentum dependence of the 
dominant peaks in the dynamical spin structure factor, 
we choose representative directions in momentum space, 
and plot 5'" L (k, uj) along such cuts. We first consider mo- 
menta lying along the (5 2) direction, shown in the left 
of Fig. ED This cut passes through two of the major 
magnetic Bragg peaks, namely at |ks|/7r w 1.3, and 3.2. 
In the low-energy region, below uj/JS w 2, most spectral 
weight is located along straight lines, emerging from the 
magnetic Bragg peaks, and with spectral weight that in- 
creases for decreasing energy, characteristic of magnetic 
soft- modes. Near the magnetic Bragg peaks, kg, we thus 
observe linear dispersion relations of the magnetic exci- 
tations, u> = c|k — kg |, with an estimated spin-wave ve- 
locity c/ JS as 2.1, which is of similar order of magnitude 
than the linear spin-wave result for the square lattice 
(c/JS = y/8 « 2.83). Similar soft-modes are also ob- 
served along the cut of constant k x /ir = 2.9, shown in 



o I ' ' ' ' 1 o I ' ' ' ' 

012345012345 
Ifl/lt along (5 2) «y 1 11 along (2.9 7t , k y ) 

FIG. 13: Intensity plot of the dynamical spin structure fac- 
tor, S(k, uj), along the (5, 2) momentum space direction (left) 
and along the line (2.97T, fc y ) (right) for the S = 1/2 Heisen- 
berg antiferromagnet on the 1393 sites approximant of the 
octagonal tiling. 



the right of Fig. 1131 including a further magnetic Bragg 
peak at k# « (2.9-7T, 4.87r). The emergence of these lin- 
ear low-energy dispersion relations is consistent with the 
linear low-energy DOS found in Sec. IIV Bl and further- 
more explains the ring-like structures in S ,_L (k), as seen 
for u)/JS = 1.8 in the right of Fig. [TT] 

In spite of the absence of translational symmetry, the 
dynamical spin structure factor of the antifcrromagnetic 
quasicrystal clearly exhibits the presence of soft-modes 
near the magnetic Bragg peak positions of the quasiperi- 
odic crystal. We expect such a generic feature of mag- 
netic long-range order to be present also in other mag- 
netically ordered quasicrystals. 

In contrast, at high frequencies, uj/JS > 5.5, we do 
not observe any significant dispersion of the spectral 
weight distribution. Instead, we find two flat bands 
of only slightly modulated spectral weight located near 
uj/JS w 6.4, and 7.3 in Fig. 1141 which correspond to the 
two isolated pockets of the spin- wave DOS in Fig. |SJ We 
consider this observation as further indication for the lim- 
ited spatial extent of the corresponding eigenstates, con- 
cluded in Sec. IIV Bl from the finite size scaling behavior 
of the inverse participation ratio. 

In the intermediate energy regime, between uj/ JS ~ 2, 
and uj/JS rs 5.5, the distribution of spectral weight 
is more complex, and can be accounted for by band- 
like segments, which recur at different energies and with 
varying bandwidths, as shown in Fig. 1141 At points of 
increased spectral weight, such as for uj/JS = 3 near 
|k|/7r a 1.8, and 4.2, in the left of Fig. ^] the correspond- 
ing gapped modes furthermore show a linear dispersion 
relation. In addition, we find bifurcations emerging as 
branches of these band-segments extending towards low- 
energies. This is seen for example in the right part of 
Fig. for momenta |k| < ir. Such self-similar struc- 
tures might have been expected to dominate the dynam- 
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FIG. 14: Intensity plot of the dynamical spin structure fac- 
tor, S(k, ui), along the (1,0) momentum space direction (left) 
and the (3, 2) direction (right) for the S = 1/2 Heisenberg an- 
tiferromagnet on the 1393 sites approximant of the octagonal 
tiling. 



ical spin structure factor, due to the geometric properties 
of the octagonal tiling, reflecting its inflation symmetry. 
Nevertheless, we find them well separated in energy from 
more conventional low-energy features, that reflect the 
magnetic order in this system. 



D. Local Dynamical Spin Susceptibility 

While the dynamical spin structure factor thus exhibits 
the peculiar nature of the excitations in the self-similar 
quasiperiodic system, the different local environments of 
the magnetic moments in the quasicrystal are accessible 
from the local dynamical spin susceptibility, the imagi- 
nary part of which at each lattice site i is given by 

xLalihUl) = S ± (i,i,C0) = (l^fel 2 + l^l 2 ) S(iO-LO k ), 

k 

(39) 

within linear spin-wave theory. The local dynamical spin 
susceptibility is accessible in nuclear magnetic resonance 
experiments, in the form of Knight-shifts at nuclear sites 
in the vicinity of the magnetic sites. Here, we study 
the properties of x[' OC ai f° r tnc Heisenberg model on the 
octagonal tiling, as an example of a quasiperiodic lattice 
structure. 

In Fig. 1151 Xiocai * s shown for the largest approximant, 
averaged separately for sites of different coordination. 
We observe a broad spread in the signal, with character- 
istic energies that increase linearly from w / JS ~ 2.5 for 
threefold sites to oj/JS w 7.3 for eightfold sites. Further- 
more, the energy-range over which there is a large signal 
narrows for sites of increasing coordination, reflecting a 
similar trend in the spread of the local staggered mo- 
ments (c.f. Fig. |3J). The inequivalent local environments 
of the various sites are responsible for this extended range 



FIG. 15: Imaginary part of the local dynamical spin sus- 
ceptibility, x'i'ocaiih w )i f° r t ne 5=1/2 Heisenberg antiferro- 
magnet on the octagonal tiling, averaged separately over sites 
with coordination numbers z = 3 to 8. The left inset shows 
Xiocaiih f° r the fivefold sites, averaged separately over sites 
with a small and large staggered magnetization, respectively. 
The right inset shows Xiocai (h 1 ^) f° r t ne eightfold sites, aver- 
aged separately over sites with different behavior under defla- 
tion transformation, grouped according to the value of z' . 



of signals. For example, the left insets of Fig. Il~5l exhibits 
that the two types of fivefold sites show a different fre- 
quency dependence of the dynamical spin susceptibility. 
Namely, sites with a smaller moment have signals inside 
a narrow region near to/JS ~ 4.5, whereas sites with 
a larger moment produce signals over a more extended 
region, ranging from u/JS rs 4.2 to u/JS w 5. The 
high symmetry eightfold sites also exhibit a characteris- 
tic splitting in the local spin susceptibility, as seen in the 
right inset of Fig. Here, the individual signals are la- 
beled by the value of z' for the site from which this signal 
results. We observe from Fig. |3 that the widths of the 
signals narrow towards u/JS ~ 7.3 for increasing values 
of z' . The local dynamical spin susceptibility thus reflects 
the hierarchical structure of the local moment distribu- 
tion of the eightfold sites. We expect such features to be 
generic properties of quasiperiodic magnets, eventually 
seen in nuclear magnetic resonance experiments on real 
quasiperiodic magnetic systems. 



V. CONCLUSIONS 

We studied the antiferromagnetic spin- 1/2 Heisen- 
berg model on the octagonal tiling, a two-dimensional 
quasiperiodic lattice structure, using linear spin-wave 
theory in a real space formulation. This approach 
was found to quantitatively reproduce previous quan- 
tum Monte Carlo results on static magnetic ground state 
properties of this system. The spin- wave excitation spec- 
trum was found to consist of magnon-like low-energy soft- 
modes with a linear dispersion relation near the magnetic 
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Bragg peaks, characteristic to long-range magnetic order. 
It will be interesting to confirm the existence of such lin- 
ear soft-modes in the octagonal tiling in future quantum 
Monte Carlo studies of the dynamical spin structure fac- 
tor. 

In addition, the dynamical spin structure factor shows 
self-similar structures and bifurcations, as well as flat 
bands at higher energies. We expect such features to be 
generic to magnetic quasicrystals, which might eventu- 
ally become observable in neutron scattering experiments 
on magnetically ordered quasicrystals. 

Within the spin- wave approach, it is possible to include 
magnetic frustration, offering the possibility of modeling 
more realistic quasiperiodic lattice structures, and their 
influence on magnetic properties. Starting in the unfrus- 
trated limit, and increasing the magnetic frustration, the 
evolution of the classical Neel state can be examined, as 
well as the potential relevance of the multifractal excita- 
tions for its breakdown in the case of strong frustration. 
Another route to magnetic disorder, which can be taken 
at least theoretically, is by means of a quasiperiodic bi- 
layer, which is expected to show a quantum phase transi- 
tion upon increasing the interlayer coupling, due to local 
singlet formation. The presence of multifractal excita- 
tions in the quasicrystal might be of possible relevance 
to quantum criticality in the transition between the Neel 
ordered state and the disordered, gaped state. We leave 
such studies for future research. 
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APPENDIX A: NUMERICAL BOGOLIUBOV 
TRANSFORMATION 

In this appendix we describe an alternative numerical 
scheme of finding the Bogoliubov transformation of the 
spin- wave Hamiltonian^i. We need to construct a matrix 
T, which diagonalizes the Hamiltonian matrix M defined 
in Eq. (JSJ, and also satisfies the constraint in Eq. (|14fl . 
We therefore need to simultaneously solve 

T f MT = O, and T f ET = T. (Al) 

The matrix T can be constructed in two steps as fol- 
lows: In a first step, an eigenvector matrix Z and the 



set of eigenvalues Xi is obtained for the non-Hermitian 
eigenvalue problem 2 ^, 

EM Z = ZA, (A2) 
where A = diag(Ai, . . . , Aj\r n ). The columns of Z are the 
right eigenvectors of EM. Hcrmitian conjugation of the 
above equation yields 

(E^EjEM = A(EZ t E), (A3) 

so that the rows of T,Z^T, form the left eigenvectors of 
SM. 

In a second step, we diagonalize the Hermitian matrix, 

L = Z^Z, (A4) 

obtaining a unitary matrix U, such that 

U^LU = diag(/i,...,ZArJ. (A5) 

Furthermore, from the eigenvalues U of L we construct 
the diagonal matrix Z^ 1 / 2 , defined as 

l^ /2 = W 1/2 - (A6) 

In case that the eigenvalues of EM are non-degenerate, 
performing this second step is trivial: In this case U is 
the E-norm of the corresponding eigenvector, and {7=1. 
We finally obtain the transformation T as: 

T = ZUr 1/2 . (A7) 

Indeed, since 

f = T f ET (A8) 

satisfies T 2 = 1 by construction, T satisfies the constraint 
Eq. (|L4|) . Furthermore, right and left eigenvectors be- 
longing to different eigenvalues are orthogonal. Thus L 
is block diagonal with blocks corresponding to degenerate 
subspaces of EM. It follows that U is block diagonal as 
well. Therefore U only mixes columns of Z belonging to 
the same eigenvalue of EM, and consequently T satisfies 

EMT = TA. (A9) 

Multiplying the above equation from the left by = 
rT*E (which follows from Eq. JSEJ), yields 

T f MT = AO. (A10) 

Hence, the matrix T also diagonalizes the Hamiltonian 
matrix M. The energies of the bosonic eigenmodes are 
given as u>i = JiST^Aj, respectively. 
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